A linear program for testing local realism 

Matthew B. Elliot 

We present a linear program tliat is capable of determining whether a set of correla- 
tions can be captured by a local realistic model. If the correlations can be described 
by such a model, the linear program outputs a joint probability distribution that 
produces the given correlations. If the correlations cannot be described under the 
assumption of local realism, the program outputs a Bell inequality violated by the 
correlations. 
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I. INTRODUCTION 



Ever since the pivotal work by Bell [hi, there has been a considerable effort to try and 
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understand the quantum-mechanical violation of local realism (as a sample, see Refs. [2, 13|, 



10. 



Bell inequalities 
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12] )■ The typical approach in most work on the subject is to search for 



13j. Experiments to test the quantum- mechanical violation of local realism 
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16|. In this paper 



have then been designed with a particular Bell inequality in mind 
we take a slightly different approach: given a set of correlations, we show how a relevant Bell 
inequality can be found. Thus, given an experimental setup, meaning the ability to create 
a particular quantum state and the ability to make a particular set of measurements, we 
illustrate how to derive a relevant Bell inequality which uses the given experimental setup 
to demonstrate a violation of local realism. 

Peres has already presented a general approach to finding all the Bell inequalities given a 
set of measurement configurations [13]. In this paper, we reproduce this result via another 
computational method involving linear programming. Specifically, we present a linear pro- 
gram which can be used to generate all the Bell inequalities for a given set of measurement 



configurations. Linear programs constitute a c 
lems for which many efficient algorithms exist 



ass of widely applicable optimization prob- 



171]. Moreover, there exists a rich theory of 
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linear programming which we can exploit. 

Our linear programming approach to problems involving local realism goes beyond finding 
all the Bell inequalities. Given an experimental setup, our linear program is capable of 
honing in on a single relevant Bell inequality. The linear program does not simply list all 
the Bell inequalities, of which there may be a large number, and present one that is violated. 
Instead, it directly finds a Bell inequality that can be used to demonstrate a violation of 
local realism. Thus, rather than tailor an experiment to the available Bell inequalities, one 
can use our linear program to find a Bell inequality tailored to an experiment that can be 
carried out. This allows the possibility of using a more robust and experimentally feasible 



setup for w^ 



Refs. 



lich a relevant Bell inequality can then be found. 
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191] also employ a linear program to find Bell inequalities. One major difference 
in our approach is in the event that a particular experimental setup is incapable of demon- 
strating a violation of local realism. In that case, no Bell inequality exists which is violated 
by the observed correlations, and our linear programming technique generates a local realis- 
tic description of all the observed correlations in the experiment. This is useful, for instance, 
in challenging claims of relevant Bell inequalities [2^. We now move to demonstrate how 
all of this is accomplished. 

Suppose, for simplicity, that we have two separated systems on which measurements of 
Aj, 1 ^ J ^ ^A, and Bk, 1 < k < itlb, are made. The generalization to more than two 
systems is straightforward. After the measurements are made, we can then compare the 
results and determine the experimentally observed correlations, (AjBk). If we allow the 
possibility of no measurement, meaning that Aj = I or Bk = I , then these average values 
also include single system average values. If no measurement is made on either system, 
meaning Aj = I and Bf^ = I, then we obtain the normalization condition. 

The existence of a local realistic description of the correlations depends on the existence 
of a joint probability distribution, 

P{Ai ^ ai, . . . , Am^ am^l Bi ^ bi, . . . , B^b ^ brr^g) = Pi{Aj aj}; {Bj bj}), (1) 

that gives the probability that a measurement of Aj yields the result aj and a measurement 
of Bk yields the result bk- This probability distribution must be normalized, 

J2 E Pi{A,--a,};{Bk^bk}) = l, (2) 
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and non negative, 

ya,,h Pi{A,^a,};{Bk^bk})>0. (3) 
It should also reproduce the observed correlations, 

ai,...,am^ bi,...,b,ng 

The existence of this joint probability distribution is equivalent to the existence of a local 
realistic description of the observed correlations. 

II. LINEAR PROGRAMMING APPROACH TO LOCAL REALISM 

In order to see the applicability of linear programming to questions about local realism, 
we consider a column vector, P, whose entries are the probabilities in Eq. ([T]). In the common 
scenario of the measurements having two outcomes, this vector has dimension 2"^-*+"^^. In 
terms of this vector notation for the joint measurement probabilities, the normalization and 
non-negativity constraints become, respectively, 

l^P = 1 and P > 0, (5) 

where 1 and are vectors of appropriate dimension with all entries equal to 1 or 0, and we 
use the linear programming convention that A> B when each entry of A is greater than or 
equal to the corresponding entry of B. 

According to Eq. (jl]) , any average value can be written as a linear combination of entries 
in P. Let us write the average values that can be observed in some experiment, of which 
there will be (m^ + l){mB + 1) if we include no measurement, as entries in a vector, C. 
Then this implies that a local realistic description of these correlations exists if and only 
if, for the appropriate matrix M, there exists a vector P with non-negative entries such 
that MP = C. In the case that measurements have two outcomes ±1, M will be an 
{rriA + + 1) X 2™-*+™^ matrix whose entries are ±1. In general, M is completely 

determined by the number of measurements and the number of outcomes, as well as their 
allowed values. Thus the problem is, given M and C, to solve 



MP = C s.t. P > 0. 



(6) 
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Notice we do not include the normalization constraint separately because it is contained in 
the equation MP = C. In fact, we assume that the first row of M has all I's and that 
the first entry in C is 1, as this takes care of normalization. Problems like these have been 
studied extensively and are called linear programs. 

A. Linear programming 

— * — * 

The most general form of a linear program is, given L, M, and C, to solve 

minimize L^P s.t. MP > (5 and P > 0. (7) 

Just about any linear optimization problem involving linear constraints can be put into this 
form. We will see later how our problem can be put into this form, but first there is an 
important feature of linear programming needing discussion. For every linear program there 
is a related dual linear program, 

maximize C^Q s.t. M^Q < L and Q > 0. (8) 

Notice that P appears nowhere in the dual program. This program seeks a completely dif- 
ferent quantity Q that generally has a different dimension than P. Yet there is an important 
relationship between these problems. The strong duality theorem for linear programming 
states that the minimum value of L^P in Eq. ([7]) is exactly equal to the maximum value 
of C^Q in Eq. (IHl) [17]. We will see that duality in linear programming reveals a duality 
between probability distributions and Bell inequalities. 

We now make a few comments regarding the theory of linear programs. First, the com- 
plexity of the optimization problems in Eqs. ([7]) and ([8]) is polynomial in the dimensions of 
L, M, and C. Also, the optimal solutions found are not unique, as there may be several 
vectors achieving the optimal value for a given linear program. 

B. A linear program for probability distributions 

Returning to Eq. ([6]), we try to rewrite our problem as a linear program, 

minimize 1 ^(MP - C) s.t. MP > (5 and P > 0. (9) 
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— # — * 

A solution to Eq. ([6]) exists if and only if this minimum is 0. The reason is that MP > C 
assures that every entry of MP — C is non negative. Thus the sum of all the entries in 
MP — C is if and only if each entry of MP — C is zero, i.e MP = C. Now we easily turn 
this into a linear program, 

minimize (M^l fP s.t. MP >C and P >6. (10) 

If the minimum value is 1 ^C, then there exists a probability distribution, P, producing the 
given measurement correlations, C. Otherwise, no such probability distribution exists. 

This linear program gives us a method for determining whether a set of correlations 
violates local realism. In addition, if the correlations can be reproduced in a local and 
realistic way, then the output, P, of the program constitutes such a description. 

C. A linear program for Bell inequalities 

There is a deep connection between probability distributions and Bell inequalities that 
arises from duality in linear programming. The dual program of Eq. (ITOl) . which produces 
Bell inequalities, is 

maximize C^Q s.t. M^Q < M^l and Q>0. (H) 

From the strong duality theorem we know that there exists a local realistic description of 
the correlations in C if and only if this maximum value is 1 '^C. Since Q = 1 satisfies the 
constraints and achieves this value, we know that generally the maximum value must exceed 

From this we deduce that there exists a local realistic description for a set of correlations 
if and only if the linear program, 

maximize C^{Q - 1) s.t. M^(Q - 1 ) < and (Q - 1) > -1 (12) 

outputs 0. This gives us the set of Bell inequalities C^{Q — 1) < for all vectors Q — 1 
satisfying the constraints. This is a complete set of Bell inequalities in that they are both 
necessary and sufficient for a local realistic description to exist. They all have the form of a 
linear combination of correlations being less than or equal to 0. 

It may appear that Eq. (fT2|) constitutes an infinite number of Bell inequalities. However, 
notice that the vectors satisfying the constraints form a convex set with a finite number of 
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extreme points. These extreme points, by themselves, then constitute a finite set of Bell 
inequalities which are necessary and sufficient for the given correlations to be described in a 
local realistic framework. Thus we do not have an infinite number of Bell inequalities, but 
only a number equal to the number of extremal vectors in the constraint set. Moreover, this 
complete set of Bell inequalities is independent of C, the correlations, and can be determined 
given only the matrix M, which is in turn only a function of the number of measurements 
as well as the possible measurement outcomes. 

If we run the dual program in Eq. (ITT!) and find a vector Q that satisfies the constraints, 
but such that C^Q > PC, then we have a Bell inequality, C^iQ — 1) < 0, that is violated 
by the given correlations in C. At this point we have completed what we set out to do: given 
a set of correlations, we have a linear program that either finds a local realistic description 
of them or produces a Bell inequality that they violate. 

As a final exercise, we show how these Bell inequalities can be simplified to yield a 
standard form for Bell inequalities. The first thing is to recognize that the constraint {Q — 
1) > — 1 is not needed. For, suppose that we have a vector Q — 1 satisfying M^{Q — 1 ) < , 
but such that {Q — 1) > —ql for some g > 1. Then the vector {Q — l)/q satisfies all 
constraints, and thus we get the Bell inequality C^{Q — 1)/^ < <(=^ C^{Q — 1) < 0. So in 
fact we get a Bell inequality for all vectors Q — 1 satisfying M^{Q — 1) < 0, regardless of 
whether they satisfy the second constraint. 

A second simplification can also be used to deal with the constraint M'^(Q — 1 ) < 0. 
Since we assume the first row of M is all I's, corresponding to the normalization of P, we 
can break up Q — 1 as Q — 1 = ^— go 5 where q has one less dimension than Q. If we also 

break up (7 = ^ 1 c j , then the Bell inequalities are —qo + c'^q < for all —qol + M'^q < 0, 
where M is M without the first row. To have a complete set of Bell inequalities it suffices to 
pick the smallest go and so we have go = (M^g)max, where the maximum denotes the largest 
entry of the vector, giving us all the Bell inequalities, 

c^q< {M^q)m.^ for all q. (13) 

This is often how Bell inequalities are stated; all the Bell inequalities are obtained by con- 
sidering an arbitrary linear combination of correlations and then saying that this must be 
less than the maximum possible value of those correlations under extremal local realistic 
predictions. 



7 



III. CONCLUSION 

What we have shown here is a computational method that determines whether a set of 
correlations can be described within a local realistic framework. If the set is so describable, 
then the linear program gives a probability distribution resulting in those correlations. If 
the set is not describable by a local realistic model, then the linear program is capable of 
generating a Bell inequality that the correlations violate. 

As a practical matter, this program can be used to generate a relevant Bell inequality 
for a particular Bell test experiment. In other words, it takes some of the difficulty away 
from designing an experiment to test the quantum violation of local realism and puts that 
difficulty into a program that yields a relevant Bell inequality to study. This is particularly 
useful in the multipartite scenario where relevant Bell inequalities are less plentiful. In terms 
of theory, linear programming gives a way to generate all the Bell inequalities and through 
the principle of duality in linear programming we see that probability distributions and Bell 
inequalities are also dual to each other. 

We note that the linear program presented here scales exponentially in the number of 
systems and the number of measurements made on the systems. This, of course, is not 



surprising because this problem is known to be NP hard [13|. For a reasonable number 
of systems and measurements, however, such as could be managed in an experiment, this 
program should still be of use. 

As a final observation, since the source of the correlations, C, is irrelevant, our linear pro- 
gram is also applicable to correlations that go beyond measurements on quantum states. In 



fact, one can consider the framework of generalized probabilistic theories |2ll], which encom- 



22|, 



passes not only quantum correlations but also the stronger correlations of PR boxes 
and determine whether a set of given correlations can be reproduced in a local realistic 
framework. 
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